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ABSTRACT 

The presence of a well-defined and narrow dust lane in an edge-on spiral galaxy is the observa- 
tional signature of a thin and dense molecular disk, in which gravitational collapse has overcome 
turbulence. Using a sample of galaxies out to z~l extracted from the COSMOS survey, we 
identify the fraction of massive (Ly) disks that display a dust lane. Our goal is to explore the 
evolution in the stability of the molecular ISM disks in spiral galaxies over a cosmic timescale. We 
check the reliability of our morphological classifications against changes in restframe wavelength, 
resolution, and cosmic dimming with (artificially redshifted) images of local galaxies from the 
Sloan Digital Sky Survey. 

We find that the fraction of Ly disks with dust lanes in COSMOS is consistent with the 
local fraction 80%) out to 2; ^ 0.7. At z=0.8, the dust lane fraction is only slightly lower. A 
somewhat lower dust lane fraction in starbursting galaxies tentatively supports the notion that a 
high specific star formation rate can efficiently destroy or inhibit a dense molecular disk. A small 
subsample of higher redshift COSMOS galaxies display low internal reddening (E[B-V]), as well 
as a low incidence of dust lanes. These may be disks in which the growth of the dusty ISM disk 
lags behind that of the stellar disk. We note that at z=0.8, the most massive galaxies display 
a lower dust lane fraction than lower mass galaxies. A small contribution of recent mergers or 
starbursts to this most massive population may be responsible. 

The fact that the fraction of galaxies with dust lanes in COSMOS is consistent with little 
or no evolution implies that models to explain the Spectral Energy Distribution or the host 
galaxy dust extinction of supernovae based on local galaxies are still applicable to higher redshift 
spirals. It also suggests that dust lanes are long lived phenomena or can be reformed over very 
short time-scales. 
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to the vertical stability of the ISM. Less massive 
spiral disks show a much more vertically fractured 
dust morphology, as turbulence, driven in part 
by stellar winds and super novae, wins out over 
the gravitational collapse towards the mid-plane. 
Specifically, D04 found that bulgeless undisturbed 
spirals with Vrot > 120 km/s consistently display 
dust lanes, while less massive disks do not. D04 
selected their galaxy sample to be "pure" disks, 
i.e., without a spheroidal stellar bulge. When 
one considers galaxies with bulges a s well, dust 



disturbances that prevent the stellar orbits from 
forming an organized bar structure. A dust lane, 
on the other hand, is an indication that the ISM 
disk is in a vertical equilibrium, no longer domi- 
nated by turbulence. Both phenomena, therefore 
indicate that the disk had some time to settle. 

As for the evolution of the bar fraction with 
redshift, s ome studies find it to be constant 
since z~l (jjogee et al. 2004) . However, there is 



lanes become common at 150 km/s (lObric et al. 
2006[ ). Dust lane prevalence in massive edge-on 



disks can be seen in the reddening of their stellar 



populations (e.g., ISeth et al.l 120051 ). or the rela- 



tive sym nietry of star-forniation tracers (24 jam 



and Ha, iKamphuis et al.l l2007l). as well as di 
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j2005l:lYoachim & DalcantonI 


20061: IWatson et al. 
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The results from 



these studies are consistent with models of the 
Spect ral Energy Distribution of these systems 
(e.g., IPodcscu et al."2000'; 'Popescu Tuffs"2005: 



Pierini et al. 2004; Baes et al. 2003, 2010; Bia nchi et al 
20001 : lBianchill2007L 120081 : iJonsson et al.ll201o] ) 



Many of the morphological features of spiral 
disks can be expected to change in the course of 
their (secular) evolution, including the scale length 
(iRavindranath et al.ll2004lTruiillo & Asuerri'2004); 
Kanw ar et al. 2008 ), surfac e brigh tn ess (Barden et al.t 



20051: ICameron fc Driverl \200li IScarlata et al 



2007 ) , and truncation radi us dTrujillo 



20051 : lAzzollini et al.l I2008D of the 



Pohlen 
stellar disk 



as w e ll as the bulge- to- dis k ratio (j Allen et al 



20061: IP annella et al.' "2009'), and bar fraction 



some controversy whether this result may be bi- 
ased by the inclusion of less-massive (and hence 
less dynamically evolved) disks in the local sam- 
ples (jBarazia^^aDlio^ By 
studying the dust lane fraction instead, we aim 
to explore an alternative diagnostic to settle this 
question. 

From first principles, equally plausible argu- 
ments can be made to expect a higher, a lower, 
or a constant dust lane fraction at earlier epochs. 
On the one hand, one might expect fewer dust 
lanes in massive spiral disks at high redshift be- 
cause (a) specific star-formation, the driving force 
of ISM turbulence, was a factor of three higher 
at higher redshift for a given stellar mass galaxy 
(jNoeske et al.l bOOTal lbllTr which would then lead 
to thicker ISM disks, and (b) molecular and dusty 
disks are still growing and may not yet match 
their stellar counterpart in size, so that they are 
not easily recognized against the stellar back- 
ground. On the other hand, there are compa- 
rable arguments as to why the dust lane frac- 
tion should stay constant or even decline with 
time. First, massive spirals should already be 
evolved enough at z~l to have their dusty ISM 



(iJogee e t al. 2004; Barazza et al. 2008; She th et al, 
120081 ). All these morphological features of the stel- 
lar disk have been explored using deep imaging 
surveys with the Hubble Space Telescope {HST). 

The results from these studies demonstrate 
that, over time, massive disks grow moderately 
in size (both the scale-length and the truncation 
radius), dim substantially in surface brightness (~ 
1 mag), and acquire bulges. The dynamical matu- 
rity of a spiral disk at high redshift - i.e., whether 
or not it has settled in an equilibrium- can be 
explored via two morphological phenomena: the 
presence of a stellar bar or that of a dust lane. The 
presence of a stellar bar is an indication that the 
stellar disk is no longer subjected to dynamical 



collapsed in to a thin di sk (e.g., iBournaud et al. 



[2009: Martig et al.llioTil ). Secondly, the fragmen- 
tation of a dust lane due to turbulence predomi- 
nantly occurs in low-mass haloes (D04), implying 
that the dust lane fraction for massive Ly galax- 
ies remains constant. And finally, a dense molecu- 
lar phase is predicted to dominate higher redshift 
galaxies ( Obreschkow fc Rawlingsl l2009bl ). ensur- 



ing the presence of a dense ISM from which dust 
lanes can form early on. 

This paper seeks to explore the possible evolu- 



-"^The cosmic star-formation density famously drops by an 
order of magnitude over the last 7-8 Gyr ago (z~l(see 

iHopkins et al.l l2008l , for the latest compilation of the 
"Madau" plot). 
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tion in the dust lane fraction. In §§[2] and O we 
discuss the data and dust lane classification. In 
§§[11 and [3 we present the derived dust lane statis- 
tics and a first-order check for observational biases. 
We analyze, discuss and summarize our results in 
§§[6l[lland[8l 

2. Data 

We select edge-on Ly galaxies from two public 
catalogs. Our primary galaxy samp le is extracted 



from the HST COSMOS survey (jScoville et al 
2QQ7al ). For comparison, and as a local refer- 
ence sample, we use the seventh data release 
from the Sloan Digita l Sky Survey (SDSS-DR7, 
2QQ9[ ). Our sample selection 
those for face-on 



[Ab azaiian e t al, 
criteria are similar to 



sive spirals used by LSheth et all (|2008h but with 
axis ratios set to select fiat galaxies (similar to 
Karachentsev et al...l993i : .Kautsch.2QQ9.) . 
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Fig. 1. — The major (A) over minor {B) axis ra- 
tio derived by Source Extractor (SE) and by a 2D 
fit (GIM2D) for the ACS observations of COS- 
MOS galaxies. The dotted line marks our selection 
criterion of A/B > 4, and the dashed line indi- 
cates agreement between both methods. The dif- 
fere nt symbols ind i cate t he classifications derived 
by (|Scarlata et al.lEoQTl ). they demonstrate that 
the vast majority of systems are indeed highly in- 
clined galaxies. 



2.1. COSMOS 

For our sample of highly inclined galax ies in 
the COSMOS survey (|Scoville et al.ll2QQ7allbl ). we 
combine information from the ACS catalog of 



the COSMOS field (ICapak et al.l 120081). th e pho- 
tometric redshift catalog (|Mobasher et al 
and the Zurich Estim ator o f Spiral Type catalog 
(ZEST, IScarlata et IDl2007h . Our main require- 
ments are that the axis ratio listed in lCapak et al.l 
(2008) is A/B > 4, and rest-frame V-band abso- 
lute magnitude My < — 21.7 in the Mobasher et 
al. catalog. This is a fixed cut in absolute luminos- 
ity to ensure a complete sample at every redshift. 
The magnitude limit is set brighter than for the 
SDSS, because an Ly galaxy fades by about one 
magnitude since z = 0.9 (Capak 2004). 

We check the reliability of the axis ratio from 
Source Extractor catalog of the ACS fields to 
select bona-fide edge-on spiral galaxies in two 
ways: first. Figure [T] compares the axis ratio es- 
timated by Source Extractor ( Bertin fc ArnoutsI 
19961 : iHolwerd a 2 005) to that of a more detailed 
fit using GIM2D (|Simardlll998l ). In general, both 
methods agree reasonably well, confirming that 
these are bona-fide fiat galaxies. Secondly, we in- 
dicate in Figure [T] the ZEST morphological clas- 
sifications which classify almost all of our galax- 
ies as highly-inclined or edge-on spiral galaxies. 
Only two "face-on" systems were originally se- 
lected, these turn out to be interacting systems 
and were excluded from the sample. 

The COSMOS survey is complete for Ly galax- 
ies out to z ~ 1, and at this redshift covers a sim- 
ilar volume as the whole SDSS survey. Therefore, 
the statistics of massive edge-on spirals in COS- 
MOS should be reliable, at least for the higher 
redshift bins (see Table [2j). A redshift of z ^ 0.9 
is the practical limit to our search, because at this 
distance, the restframe wavelength of COSMOS 
F814W filter corresponds to SDSS-^ filter which 
is the bluest usable SDSS passband (see § 12. 2p . 

We point out that an HST pixel corresponds 
to 0.4 kpc at this distance, which is comparable to 
the typical scale heig:ht of dust lanes (^ 0.5 kpc, see 
Seth et al.ir2005l : lKamphuis et al.ll2007f ). so that we 
can readily identify dust lanes in the COSMOS 
images. 
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2.2. SDSS 

SDSS galaxies were selected by oblateness 
{A/B > 4), major axis {A > 30'', to ensure 
the selection of well-resolved objects), and for 
My < — 20.7, using a combination of g-r color 
and ^-magnitude selectiorjl. The corresponding 
SQL quer}|f| yields 194 galaxies. 

We used this SDSS sample to simulate the ap- 
pearance of local galaxies at high redshift. For 
this, we produced a "redshifted" SDSS catalog 
(z-SDSS hereafter) by rebinning the u^g^r and 
i images of all galaxies, cosmologically dimming 
them, and adding appro priate sky nois j j, follow- 



ing the prescription from iGiavalisco et al.l (|1996f ). 
We used the SDSS Atlas images, i.e., those with 
sky-pixels and other objects in the field set to 
blank and added the appropriate ACS sky noise. 
The SDSS-2; images were not redshifted as they 
do not overlap with the COSMOS F8I4W fdtei at 
any redshift other than zer . We use this artifi- 
cial z-SDSS catalog to check the COSMOS results 
for observational biases, as discussed in more de- 
tail in gSl 

3. Identification of Dust lanes 

One of us (BWH) identified dust lanes by eye 
in single filter images, after randomizing the orig- 
inal five filters of the SDSS, the four filters of 
redshifted SDSS images (see g5j), and COSMOS 
images before presenting them for classification. 
This method provides an objective assessment of 
whether or not a galaxy has a dust lane, because 
the classifier is not infiuenced by a previous dust 
lane classification in a bluer filter in SDSS or lower 
redshift in the redshifted SDSS images. Using a 



■^Conversion to Johnson V-magnitude from 
SDSS was done using the tools available at 
[http: //www. aerith.net/astro/color_conversion.html , 
and the SDSS f Cosmo Abs Mag function, using t he photo- 
metric system described in lFukugita et all (|l99Q V 

3 G.isoA.g/G.isoB.g > 4, G.isoA.g > 30/04 ^ND S.z > 
AND dbo.fCosmoAbsMag((G.r+1.48*(G.g-G.r)+0.17), 
S.z) < -20.7 . 

^The SDSS input pixelscale was set to 0(^9 in order to ac- 
count for the different resolutions of the SDSS and ACS 
images; the resolution difference is a factor 20 while the 
difference in pixel scales is 8. 

^HST/ACS's filters include the exact sdss-z: F775W but 
COSMOS observing strategy opted for the widest red filter 
available wide-/, F8I4W 



combination of zoom, rotation and change in con- 
trast, each galaxy was visually classified for the 
presence of a dark band in the plane of the disk. 
Identification is facilitated by the fact that a dust 
lane is a galaxy-wide phenomenon, but becomes 
rapidly more difficult as the inclination strays from 
edge-on (90°) inclination. 

Post facto, we checked the internal consistency 
of the classification as a function of wavelength 
for the native and redshifted SDSS sample; i.e., 
whether the number of identified dust lanes agrees 
with those for SDSS-^ filter which serves as the lo- 
cal reference. Three filters (SDSS-^; r and i) are 
internally very consistent (agreement more than 
80%, Table [1]). The classifications in the redshifted 
SDSS g, r and i images (z-SDSS) are also very 
consistent with the zero-redshift SDSS-^ classifi- 
cations, as well as internally consistent with each 
other (better than 80%). Thus, our classifications 
are reliable to 80% or better for the SDSS-^^ rand 
i filters or over the redshift range 0.1 to 0.8. 

While we did classify the sdss-u images, both 
original and redshifted images (Table [1]), the re- 
sulting dust lane fractions are not usable as these 
images have too poor signal-to-noise to reliably 
identify dust lanes and, once redshifted, a ~lkpc 
wide dust lane is no longer resolved. 

Alternative ways to identify dust lanes, a ver- 
tical color profile or a color composite image, 
are currently not possible with the COSMOS 
data as only a single filter (F814W) is avail- 
able for the HST observations. However, in the 
future, the CANDELS (The Cosmic Assembly 
Near-infrared Deep Ex t rag:alactic Legacy S ur vey 
iKoekemoer et al.l I2OIII : ICrodn et alTl201l[ ). wih 



provide near-infrared and optical images for a 
subsection of this and several other fields and a 
color-based analysis could be done with uniform 
data and better statistics. 

4. Initial Results 

Figure [2] shows, as a function of redshift, the 
fraction of the COSMOS sample that shows a dust 
lane. The error est imate is based on equation 
26 in iGehrelsl (|T986[ ) with a confidence level of 
95%, for small numbers in a redshift bin {N < 
30), and for larger numbers, we use simply a = 
1/ (/(I — f)/N), where / is the dust lane fraction 
and TV is the total number of galaxies in the red- 
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Fig. 2. — The observed dust lane fraction in lu- 
minous galaxies (brighter than Ly) for the SDSS 
(triangles, green: SDSS-^, blue: SDSS-r, and 
red: SDSS-i, respectively) and COSMOS galaxies 
(black circles). These fractions are listed in Ta- 
ble [2l These data are yet unmatched in rest-frame 
wavelength or resolution, and cannot be used to 
correctly infer evolution in the dust lane fraction. 



20081 ) . 



shift bin (see also iDesai et al.l I2QQ7I : iSheth et al 



The error is largest in the lowest redshift bin 
because COSMOS probes a smaller volume (and 
hence lower number of objects) at low redshift 
(see Table [2]). We also show the dust lane frac- 
tion in the SDSS-^ filter as the local reference. 
This fraction is not as high as the one found by 
D04 for flat, bulgeless, massive galaxies (~ 100%), 
but is consistent with t hat for a mix of Hubble 
types tobri c et al.ll2006[ ). The SDSS fraction rises 
slightly if flatter (higher values for A/B) galax- 
ies are considered; i.e., if the selection favors even 
more bulgeless galaxies. Therefore, we can dis- 
count any changes in dust lane fraction due to se- 
lection effects in the oblateness criteriorjl 

Taken at face value. Figure [2] indicates a peak 



^A/B > 4 is the same in SDSS at z = and COSMOS at 
z = 0.8 as the respective resolutions are similar. However, 
at lower redshifts the COSMOS data resolves the thickness 
of disks better with consequently lower values for B. To 
select similar galaxies in SDSS, one would need a stricter 
A/B criterion. 



in the dust lane fraction around z~0.3, with a 
decrease towards both higher and lower redshift. 
However, this trend can, in principle, be caused 
by observational effects, such as the varying rest- 
frame wavelength of the COSMOS redshift bins 
(because dust extinction is more pronounced for 
shorter wavelengths), or the reduced resolution in 
images of galaxies at large distances (because a 1 
kpc thick dust lane at z=0.8 is only just resolved in 
HST images). In the following section, we use the 
redshifted SDSS data to estimate the amplitude 
of such potential systematic errors before further 
analysis and interpretation. 

5. Systematics in Dust Lane Classification 

5.1. Restframe Wavelength 

To asses the possibility that our dust lane clas- 
sification is affected by the shift in restframe wave- 
length, we classified dust lanes in each SDSS filter 
separately (see Figures [3] and [H for examples). In 
Figure m the triangles indicate the dust lane frac- 
tion in the SDSS sample as a function of wave- 
length. Each SDSS filter has been placed at the 
redshift at which the central wavelength of the 
COSMOS F814W mtei would match that of the 
SDSS filter. SDSS-^^^i^r and g correspond to rest- 
frame at z = 0,0.114,0.352, and 0.778 for the 
F^^^W filter, respectively (Table [1]). As expected, 
the dust lane fraction appears to increase towards 
bluer rest wavelengths, due to the stronger extinc- 
tion at shorter wavelengths. Therefore, this would 
cause an apparent evolution in the COSMOS sam- 
ple that is opposite to what we have found, imply- 
ing that our results are not caused by any wave- 
length dependence in the dust lane classification. 

5.2. Signal-to-Noise Ratio 

Another possible observational bias that could 
play a role are differences in the signal-to-noise 
in the images of the various SDSS filters. The 
SDSS-2; images have notably lower S/N than the 
other filters due to the combined lower transmis- 
sion of both filter and atmosphere. Therefore, the 
observed dust lane fraction is low in the SDSS- 2; 
filter due to a combination of poor signal-to-noise 
and low contrast of a dust lane in the reddest fil- 
ter. We note that the fraction for SDSS-^^ r and 
i does not change as dramatically as from SDSS-2; 
to SDSS-i Combined with the poor consistency 
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Table 1: The redshift, corresponding look-back time (1), filter (2), restframe wavelength (3), and dust lane 
fraction for the SDSS (4) and redshifted SDSS (z-SDSS, 5). The consistency checks (6,7) of the classifiers is 
expressed as the percentage of the answers (dust lane y/n?) that agrees with the original SDSS-^ answers. 
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of the SDSS-^ with SDSS-g (see 33]and Tabled]), 
we discard SDSS-2; as a reference. For the compar- 
ison with higher redshift, we show the dust lane 
fraction in SDSS-^^ r and i filters as the redshift 
zero reference in Figure O 

5.3. Spatial Resolution 

The change in rest wavelength with redshift 
may not be the only effect biasing the dust lane 
fractions. At high redshift, galaxies are less re- 
solved, which should make dust lanes more diffi- 
cult to identify. This trend would counteract the 
bias due to bluer rest wavelengths. A minor effect 
could be that at lower redshift s, the slight smooth- 
ing of the image might aid in the identification of 
a galaxy- wide phenomenon such as a dust lane. 

To asses the total impact of these potential bi- 
ases, in Figure [5l we also plot the apparent dust 
lane fraction for the artificially redshifted SDSS-z 
images (squares). In addition to possible effects 
caused by changes in restframe wavelength, this 
test should also account for the change in resolu- 
tion and surface brightness. Because there is little 
or no significant change in the dust lane fraction 
between z ~ 0.1 to z ~ 0.8, we conclude that the 
reduced resolution and signal-to-noise at high red- 
shift at most balances the increased sensitivity to 
dust lanes at shorter wavelengths. 

5.4. Filter Bandwidth 

As a final check for the reliability of the ob- 
served evolution in the dust lane fraction, we com- 
pare the dust lane fraction in the redshifted SDSS 



images {g, r and i) to that in the COSMOS galax- 
ies whose redshifts fall within the SDSS filter's 
bandwidth. The F8I4 W filter is much wider than 
the SDSS filters when these are shifted to the red- 
shifts at which they correspond to the restframe 
wavelength. 

From the various options for defining the band- 
width, we chose to select COSMOS galaxies whose 
photometric redshifts (photo-z) falls inside the 
bandwidth of the SDSS filter (after it was red- 
shifted from its central wavelength to the central 
wavelength of the F8I4 W filter. Figure [6] shows 
the results after resampling the COSMOS classi- 
fications in this way; black symbols for those in- 
tervals for which redshifted SDSS information is 
available. The gray area is all the COSMOS sam- 
ple, including the intervening intervals. Table [2] 
lists the fraction of the COSMOS sample in the 
bins corresponding to the redshifted SDSS filters 
and dust lane fraction in the intervening intervals. 

A constant fraction of 80% dust lanes is consis- 
tent with almost the entire redshift range probed 
by COSMOS. Only close to z ^ 0.75 does the 
COSMOS fraction deviate by more than a stan- 
dard deviation. The constant fraction is consistent 
with the redshifted SDSS images (within 95% of 
each point). Thus, any evolution in the dust lane 
fraction is with only with 95% (la) confidence but 
the slight rise with redshift appears unlikely due 
to the previously discussed observational biases. 

One remaining observational biases is that the 
different redshifts in COSMOS predominantly 
probe different environments, i.e., field versus 
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Fig. 3.— Examples of the SDSS (top row), redshifted SDSS (middle row) and COSMOS (bottom) edge-on 
galaxies with dust lanes. The COSMOS examples are chosen to be as close to the SDSS redshift bin centre 
as possible. Because no COSMOS objects were selected above ^ = 1, there is no comparison object for the 
SDSS- 2/ band. 



group or cluster. If the dust lane fraction de- 
pends on environment but not secular evolution, 
this could mimick the redshift evolution we see in 
Figure [6l 



6. Analysis 

The COSMO S phot ometric redshift catalog of 
Mobasher et al. uses fits to the Spectral 

Energy Distribution (SED) to provide several use- 
ful parameters for the COSMOS galaxy sample, 
in particular the photometric type {Tphot)^ inter- 
nal reddening {E[B — V])^ and stellar mass (M*). 
Here, we explore possible dependencies of the dust 
lane fraction on these parameters. We apply the 
same bins as Figure [6] and Table [21 

The Tphot parameter is a photometric Hubble 
type, separating red, early-type galaxies {Tphot 
= 1) from the bluest, late- type g:a laxies (T^hot = 
6). The first four Tphot types are IColeman et al 



( 198Qr) templ a tes, and the last two are from 



Kinnev et al 



these are starbursts with 
little internal extinction. In Figure [71 we compare 
the dust lane fraction of quiescent {Tphot < 4) to 
that of star- forming {Tphot ^ 4) galaxies. The 
latter have sufficiently accurate statistics only in 
the three highest redshift bins. 

We note that the SED fit allows only for a 
small amount of re ddening of the stellar t emplate 
{E{B -V) < Q.25. lMobasher et aLllioQ?! ). and in 
the case of edge-on spirals, one expect much higher 
levels of extinction. However, for a distant spiral 
with a clear dust lane, the SED fit is essentially 
made to the mostly undimmed extra-planar stel- 
lar population and the mid-plane population is too 
dim to be accounted for. Therefore, some of the 
quiescent SED galaxies may, in fact, be hosting 
some star-forming regions deeply embedded in the 
dusty structures. 

The more quiescent SED galaxies {Tphot < 4), 
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Fig. 4. — As Figure [3], but for galaxies without dust lanes. 



which dominate the COSMOS population, show 
the gradual increase in dust lane fraction that is 
nevertheless consistent with a constant fraction of 
80%, similar to the full COSMOS sample and red- 
shifted SDSS images. The star-forming galaxies 
{Tphot = 4 — 6), on the other hand, have a consis- 
tently lower dust lane fraction, although the con- 
stant fraction of 80% is still within one standard 
deviation due to the much lower statistics for star- 
bursting galaxies. By combining the last three 
redshift bins in the COSMOS sample, the frac- 
tion of dustlanes in starbursting galaxies is signif- 
icantly lower than the quiescent galaxies. A lower 
dust lane fraction is consistent with the notion 
that greatly enhanced specific star-formation (a 
starburst) destroys or inhibits the formation of a 
dust lane. However, a lower level of star- formation 
can co-exist with a dust lane as the star-formation 
would inherently be deeply embedded. To give a 
definite answer, one would need a star- formation 
rate estimate based on the full SED (including 
sub- mm to estimate the embedded star-formation) 
and much improved statistics. To summarize. Fig- 



ure [71 is consistent with a constant dust lane frac- 
tion for both populations of galaxies. 

The photometric fits of iMobasher et ah ( 2007 ) 
also report the internal dust reddening of the SED 
template. Figure [8] shows the dust lane fraction 
for reddened (^[5 - V] > 0.1) and little or not- 
reddened galaxies {E[B — V] < 0.1). The majority 
of galaxies in the COSMOS sample show redden- 
ing of their SED. The small sub-sample of galax- 
ies that do not show reddening also has a lower 
incidence of dust lanes, especially at z=0.7. To 
counter the poor statistics, we combine the last 
three redshift bins in COSMOS (crosses in Figure 
[8]). The low- or no-extinction population is offset 
from the general COSMOS population with only 
30% showing a dust lane. 

This small subsample with no dust reddening 
is surprising because edge-on spirals are the ob- 
jects in which one would expect dust reddening to 
be most pronounced, with the possible exception 
of (U)LIRGS or merging systems. Yet, there is a 
small population of edge-on disks at higher red- 
shift with little reddening and no apparent dust 
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Table 2: The fraction of COSMOS and redshifted SDSS galaxies in redshift bins corresponding to the SDSS 
filter widths: (1) redshift range, (2) corresponding rest-frame SDSS filter, (3) fraction of galaxies with dust 
lanes in redshifted SDSS sample, (4) fraction of dust lanes in COSMOS galaxies, (5) 



Redshift 


0.063 - 0.089 


0.089 - 0.28 


0.28 - 0.38 


0.38 - 0.62 


0.62 - 0.95 


SDSS 


sdss-i 




sdss-r 




sdss-g 


zSDSS frac 


0.84 ± 0.03 




0.81 ± 0.03 




0.79 ± 0.03 


COSMOS frac 


1.00 ± 0.46 


0.94 ± 0.18 


0.96 ± 0.14 


0.84 ± 0.05 


0.76 ± 0.04 


Tphot 












<4 


1.00 ± 0.46 


1.00 ± 0.13 


0.96 ± 0.14 


0.88 ± 0.05 


0.78 ± 0.04 


> 4 




0.50 ± 0.44 




0.67 ± 0.30 


0.62 ± 0.23 


E(B-V) 












< 0.1 






1.00 ± 0.46 


0.50 ± 0.44 


0.27 ± 0.29 


> 0.1 


1.00 ± 0.46 


0.94 ± 0.18 


0.95 ± 0.15 


0.85 ± 0.05 


0.80 ± 0.04 


log(M/M*) 












<0.1 


1.00 ± 0.46 


1.00 ± 0.13 


0.95 ± 0.16 


0.82 ± 0.23 


0.85 ± 0.05 


> 0.1 




0.00 ± 0.45 


1.00 ± 0.33 


0.84 ± 0.06 


0.70 ± 0.05 



lane (Table [2]). If the dust lane is dispersed in these 
disks and there are dusty structures throughout 
these galaxy disks, (similar to the lower mass disks 
in D04), one would expect some reddening of the 
stellar template. We speculate that in these cases, 
the galaxy's dusty, molecular disk is much smaller 
than the stellar one, possibly after losing much of 
their ISM to the surrounding I GM through, for in- 
stance, a galactic fountain (e.g., Sturm et al.ll20lf ) 
or stripping. 

Finally, the dust lane fractions at different stel- 
lar masses are shown in Figure [9l The high and 
low mass samples are both consistent with a con- 
stant fraction of dust lanes (80%) with the excep- 
tion of the points at z=0.7. Here the higher mass 
galaxies display a lower dust lane fraction, while 
the lower mass galaxies are consistent with the 
SDSS and the constant fraction of 80%. If we av- 
erage over the last three redshift bins to improve 
statistics (crosses with no symbols in Figure [9]), 
the discrepancy between high and low-mass galax- 
ies remains. This implies a reversal of the relation 
between mass and occurrence of dust lanes seen 
in D04. However, both samples consist of massive 
galaxies as per the definition of D04 {vrot > 120 
km/s). One could speculate that the a fraction of 
the massive spirals in COSMOS at z=0.7 have re- 
cently undergone a major merger (as most massive 
spirals are predicted to have since z ^ 1). Thus, in 
some massive disks, a few of the dust lanes are now 



missing. The break in dust lane fraction observed 
by D04 occurs below a stellar mass of 0.1 M*. 
The starburst galaxies in Figure [7] could be those 
massive disks that are missing their dust lane. 

7. Discussion 

The main assumption in D04 and this paper 
is that dust traces the coldest, often molecular, 
phase of the ISM. This assum ption is based on 
both observational data (e.g ., lAllen et al.l 



19861 



Vogel et a l. 1988; Ran d 1995; Planck Collaboration et al 



2011 



Holwerda et al. submitted), and theo- 
retical arguments, which suggest that dust and 
molecular gas have slow rel ative drift velocities 
([Weingartner fc Drainell2QQl[ ). Furthermore, dust 
acts as a cataly st for molecul ar hydrogen {H2) 
production (Ca zaux fc TielensI [2004), and dust 
shields the molecular gas from dissociation by 
ultra-violet radiation. Our finding that the dust 
lane fractions seems independent of redshift out 
to z=0.7 (Figure [6]) therefore implies that the 
global state of the ISM in large spiral disks has not 
changed over the past 7 Gyr. We emphasize that 
this statement is valid only for the galaxy pop- 
ulation as a whole, and that in individual disks, 
dust lanes may well be destroyed and re-formed on 
short time-scales. Taken at face value, however, 
the constant fraction over time suggests that dust 
lanes are a long-lived phenomena. 
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A Original SDSS 
■ Redsliifted SDSS 
--- COSMOS 



sdss-g 



0.0 



0.2 



0.4 0.6 
Redsliift (z) 



COSIVIOS 
Redshifted SDSS 



0.4 0.6 

Redsliift (z) 



Fig. 5. — The fraction of galaxies with identi- 
fied dust lanes in the original SDSS images (blue 
triangle) and redshifted (dimmed, rebinned etc.) 
SDSS images as a function of redshift. The in- 
dicated redshift corresponds to the the restframe 
wavelength observed by the COSMOS F^^^W fil- 
ter. The COSMOS fractions from Figure [2] are 
shown as a dashed line with uncertainties as the 
grey shaded area. 



Fig. 6. — The fraction of dust lanes in the red- 
shifted SDSS and COSMOS, selecting only those 
COSMOS galaxies whose photo-z would place 
them in the bandwidth of the SDSS filter red- 
shifted to the central wavelength of the F8I4W 
filter. Only the redshifted SDSS-^^ r and i are 
considered. The respective widths of the shifted 
filters are marked with the horizontal errorbars. 
SDSS points are shifted by 0.005 redshift for clar- 
ity. 



The observation that approximately the same 
number of massive galaxies have thin dusty molec- 
ular disks at z=0.8 has significant implications 
for models of their emission and the supernovae 
they host. For instance, SED models of local 
massive edge-on systems are still suitable for ex- 
plaining the ch aracteristics o f these higher redshift 
syste ms 



200| 



ie^ 



Xilouris et al.j 1999: Gor don et al. 



Misselt et al 
Baes et al.' 
Holwerda et al 




Bianchil I2007L \2 



iBianchi & Xilouris' ' 2011 

Schechtman-Rook e t al.l 

2012al fl Holwerda et al. submitted). The SNIa 
lightcurve fits of Dark Energy Cosmology depend 
on a host galaxy dust extinc tion prior based on 
a flat dusty spiral disk (e.ps;.. iHatano et al.l Il998t 
Lah et al.ll2007t iHolwerdal I2OO8I ) . For those SNIa 
found in massive spiral disks (likely the majority), 
the locally derived prior would be applicable out 
to z - 0.8. 

New Herschel results suggest a strong evolu- 
tion in the population of dusty g:a laxies (e.g., 
Dve et al-llioiot IClements et al.ll2010[ ). There are 



many more sub-mm luminous galaxies at higher 
redshift then there are today, a result of both en- 
hanced dust production and higher irradiation due 
to star-formation. The Herschel results show how 
the dusty ISM's illumination rapidly evolves with 
time. This paper, in contrast, focuses on the evo- 
lution of dusty ISM's morphology in massive spiral 
disks and finds little or no evidence for evolution 
since z=0.7. Both dust illumination and geometry 
are key ingredients for models explaining the pop- 
ulation of luminous infra-red and sub-mm galaxies 
observed at high redshift. 

Because we select galaxies that ar e expected to 
passively evolve in to M* galaxies (see lCapakll2004l : 
[Sheth et al.' 2008), these systems are the massive 
haloes and disks expected to have collapse first. 
They have produced a large fracti on of their stellar 
mass at the epoch of observation ( Bell et al]l2003l : 
lllbert et al. 2010), and the mass surface density of 
all components combined (including dark matter, 
molecular and atomic gas) is in a similar vertical 
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COSMOS V,,<4 
COSMOS V„,>4 
Redshifted SDSS 



0.4 0.6 

Redshift (z) 



COSMOS E{B-V) <0.1 
COSMOS E{B-V) >0.1 
Redshifted SDSS 



0.2 



0.4 0.6 

Redshift (z) 



Fig. 7. — The dust lane fraction of edge-on 
Ly galaxies in the COSMOS survey, sorte d by 
star fo rmation type {Tphot) from Mobasher et alJ 
(|2QQ7h . Plotted are the fuh COSMOS sample 
(gray area), early- type SEDs {Tphot< 4, green cir- 
cles) and late- type SEDs {Tp^ot^ 4, blue trian- 
gles), which are extinction- free starbursts (shifted 
by z=0.01 for clarity). The numbers are listed in 
Table [2l The average for the last three redshift 
bins of COSMOS are indicated by the green and 
blue crosses without a symbol for early and late- 
type SEDs respectively. 



balance with turbulence as it is today i.e., suffi- 
cient to collapse the ISM into a thin dust lane. 
Since the specific star-formation is a factor three 
higher for the same stellar mass ( iNoeske et al. 



l2QQ7al lbl). it implies that the total ISM surface den- 
sity must be higher to compensate. The propor- 
tions of the IS M components depend on the local 
pressure (e.g., lObreschkow fc Rawling 
and we speculate that increase in ISM surface 
density and star-formation feedback both point 
to a more important molecular component at this 
epoch. 



8. Conclusions 

By identifying dust lanes in highly inclined spi- 
ral disks of LI, galaxies in the SDSS and COSMOS 
surveys, we conclude the following: 



1. Our estimate (80%) for the dust lane frac- 



Fig. 8. — The dust lane fraction of edge-on Ly 
galaxies in the COSMOS survey, sorted by in- 
ternal dust reddening of the SED as derived by 
Mobasher et al. l (|2QQ7l) . Plotted are the combined 
COSMOS sample (gray area), reddened galaxies 
{E[B — V] > 0.1, blue triangles), and galaxies with 
little or no reddening {E[B — V] < 0.1, green cir- 
cles, shifted by z=0.01 for clarity). The average 
for the last three redshift bins of COSMOS are 
indicated by the green and blue crosses without 
a symbol for low and high reddening respectively. 
At low redshift, effectively all galaxies display ev- 
idence of reddening in their SED, but there is a 
small population of galaxies at higher redshift that 
does not show any evidence for reddening or a dust 
lane. The numbers are listed in Table [2l 



tion in local galaxies, derived from SDSS-^ 
images , is consistent with previous results 
from iDalcanton et al.l (|2004l ) (Figure EJ. 



2. Any evolution in the dust lane fraction with 
redshift derived from COSMOS images can 
not be caused by the effects of varying rest- 
frame wavelength or spatial resolution for 
galaxies at higher redshift (Figure [5j). Taken 
in the same filter bandwidth, the redshifted 
SDSS and the COSMOS classifications are 
both consistent with a constant dust lane 
fraction of about 80% of galaxies since z ~ 
0.8 (Figure El). 



3. The majority of L^. galaxies in COSMOS 
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tentative. 



COSMOS log{M/M* ) <0.1 
COSMOS log{M/M*)>{).l 
Redshifted SDSS 



0.2 



0.4 0.6 

Redshift (z) 



Fig. 9. — The dust lane fraction of edge-on Ly 
galaxies in the COSMOS su rvey, so rted by stellar 
mass as derived by Mobasher et al.l (12007). Lower 
and higher stellar masses are the green circles and 
blue triangles, respectively, the latter shifted by 
z=0.01 for clarity. The average for the last three 
redshift bins of COSMOS are indicated by the 
green cross for low-mass and blue cross for high- 
mass galaxies. Both mass ranges display similar 
fractions of dust lanes out to z < 0.6, with a devi- 
ation at z=0.7 (Table E}. 



(75% or more) display a dust lane over the 
last 7 Gyr. The observed fractions in COS- 
MOS are consistent with a constant frac- 
tion of 80%, i.e., with no evolution since 
z=0.7. At z=0.8, close to the redshift limit 
of our survey, the dust lane fraction appears 
slightly lower, hinting at some evolution be- 
yond z=0.7 (Figure [6]). 

4. Star- forming galaxies (SED type T^i^ot = 4- 
6) seem to have a lower dust lane fraction 
(/ ~50%) than more quiescent SED types 
{Tphot< 4, / ~ 80%), consistent with the 
picture that a high specific star-formation 
rate injects enough energy into the ISM disk 
that the additional turbulence prevents the 
vertical collapse of the dusty molecular ISM 
into a dust lane (Figure [7]). However, we 
caution that due to the small number of star- 
forming galaxies, the associated error bars 
are rather large, so this result is somewhat 



5. A small population of edge-on spirals with 
no reddening or dust lane is evident at higher 
redshift (Table El Figure [8]). We specu- 
late that their dust disks do not extend far 
enough into the spiral disk to project a dust 
lane. 

6. At the highest redshift {z ~ 0.8), the lower 
mass galaxies {log{M /M*) < 0.1) show a 
slightly higher fraction of dust lanes than 
the high mass galaxies {log{M /M*) > 0.1, 
Figure [9|). This may be due to evolution, 
a change in environment or a higher con- 
tribution of starbursting galaxies at higher 
redshift to the massive galaxy population. 

7. Many of the Ly disks have already accumu- 
lated enough surface density 7 Gyr ago in all 
components combined to vertically collapse 
the cold, dusty ISM into a thin lane. This 
implies SED or SNIa host galaxy models us- 
ing this geometry based on local results are 
applicable to higher redshift {z < 0.8) 

Future work on the fraction of dust lanes in 
edge-on spiral disks will be possible with the CAN- 
DELS survey which will provide three color images 
of five fields. The near- infrared imaging with HST 
allows to probe further along the mass range, pos- 
sibly down to the mass where dust lanes disap- 
pear in these distant galaxies. In the Sloan Dig- 
ital Survey, we hope to use the classification by 
the GalaxyZ oo2 volunteers to improve on the re- 
sults of D04, lObric et al.l (|2QQ6h and this paper for 
the local Universe. Lastly, edge-on galaxies can be 
observed with the Atacama Large Millimetre Ar- 
ray (ALMA) to probe both the dynamics and the 
width of the molecular phase in spirals. 
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